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In this work the following lepton flavor violating r and fi decays are studied: fi~ fi^ , 

T~ e~e~e+, e~ e~ e'^ , t~ /i~7, r~ — > e~7 and fi~ e~7. We work in a con- 

strained supersymmetric scenario with universal soft breaking terms and where the MSSM 
particle content is extended by the addition of three heavy right handed Majorana neutrinos 
and their supersymmetric partners. The generation of neutrino masses is done via the seesaw 
mechanism. We analize all these lepton flavor violating decays in terms of the relevant input 
parameters, which are the usual mSUGRA parameters and the seesaw parameters. In the 
numerical analysis compatibility with the most recent experimental upper bounds on all these 
r and jj, decays, with the neutrino data, and with the present lower bounds on the supersym- 
metric particle masses are required. Two typical scenarios with degenerate and hierarchical 
heavy neutrinos are considered. We will show here that the minimal supergravity and see- 
saw parameters do get important restrictions from these r and n decays in the hierarchical 
neutrino case. 



1 Introduction 

The present strong evidence for lepton flavor changing neutrino oscillations in neutrino data im- 
plies the existence of non-zero masses for the light neutrinos, and provides the first experimental 
clue for physics beyond the Standard Model (SM). This fact does require a theoretical framework 
beyond the SM with just three massless left-handed neutrinos. Within the MSSM-seesaw con- 
text, which will be adopted here, the MSSM particle content is enlarged by three right handed 
neutrinos plus their corresponding supersymmetric (SUSY) partners, and the neutrino masses 
are generated by the seesaw mechanism. Three of the six resulting Majorana neutrinos have light 
masses, m^.. ,i = 1,2,3, and the other three have heavy masses, m^.ji = 1,2,3. These physical 
masses can be related to the Dirac mass matrix rriD, the right-handed neutrino mass matrix 
niM, and the unitary matrix Umns by diag{m,y^,miy^,my^) ~ Uj^j^g{—m£)ni^l'm^)UMNS and 
diag(m7Vi,'7iAr2,m7V3) mu, respectively. Here we have chosen an electroweak eigenstate ba- 
sis where tum and the charged lepton mass matrix are flavor diagonal, and we have assumed 
that all elements in m/5 = < H2 >, where is the neutrino Yukawa coupling matrix and 
< H2 >= f sin/3 (v = 174 GeV), are much smaller than those of tum- The two previous re- 
lations can be rewritten together in a more convenient form for the work presented here as, 

nvjy = im^^^^^'^ Rm'^^^^^^'^U^j^fg, where i? is a general complex and orthogonal 3x3 matrix, 
which will be parameterized by three complex angles 9i, i = 1,2,3. 

One of the most interesting features of the SUSY-seesaw models is the associated rich phe- 
nomenology due to the occurrence of lepton flavor violating (LFV) processes. In these models 
the LFV ratios can be large due to an important source of lepton flavor mixing in the soft- 
SUSY-breaking terms. Even in the scenarios with universal soft-SUSY-breaking parameters at 



the large energy scale associated to the SUSY breaking , the running from this scale down 
to niM induces, via the neutrino Yukawa couplings, large lepton flavor mixing in the slepton 
soft masses, and provides the so-called slepton-lepton misalignment, which in turn generates 
non-diagonal lepton flavor interactions. These interactions can induce sizable ratios in several 
LFV processes which are actually being tested experimentally with high precision and therefore 
provide a very interesting window to look for indirect SUSY signals. Among these processes, 
the LFV r and // decays are probably the most interesting ones for various reasons. On one 
hand, they get vanishing rates in the SM with massless neutrinos and highly suppressed rates 
in the SM with massive netrinos. The smallness of these rates in the non-SUSY version of the 
seesaw mechanism for neutrino mass generation is due to their suppression by inverse powers of 
the heavy scale mM- On the other hand, although these decavs have not been seen so far in the 
present experiments, there are very restrictive upper bounds^ on their possible rates which im- 
ply important restrictions on the new physics beyond the SM. These restrictions apply severely 
to the case of softly broken SUSY theories with massive neutrinos and the seesaw mechanism, 
since these give rise to higher rates, being suppressed by inverse powers of the SUSY breaking 
scale, msusY < 1 TeV, instead of inverse powers of uim- 

Here we compute the partial widths for the the LFV r and /j, decays of type Ij lij 
and Ij 3/j to one-loop order and analyze numerically the corresponding branching ratios in 
terms of the mSUGRA and seesaw parameters, namely, Mq, M1/2, tan/3, tun^ and R. To solve 

numerically the RGEs we use the Fortran code SPheno^ that we have adapted to include the 
full flavor structure of the 3x3 soft SUSY breaking mass and trilinear coupling matrices and 
of the Yukawa coupling matrices. Our final goal is to use the SUSY contributions to LFV r 
and decays as an efficient way to test the mSUGRA and seesaw parameters and we explore 
in detail the restrictions imposed from the present experimental bounds. We will find that for 
some plausible choices of the seesaw parameters, being compatible with neutrino data, there are 
indeed large excluded regions in the mSUGRA parameter space. 

For the numerical analysis, the Umns niatrix elements and the rriu- are fixed to the most 



favored values by neutrino data with y Am^^^ = 0.008 eV, y Am^^^ = 0.05 eV, 612 = 9soi = 30°, 
^23 = Gatm = 45°, ^13 = 0° and 6 = a = f3 = 0. Some results will also be presented for the 
alternative choice of small but non- vanishing ^13. We consider two plaussible scenarios, one 
with quasi-degenerate light and degenerate heavy neutrinos and with niui ~ ^-^ eV, 111^2 = 

rriy-^ + , = rriy-^ H — gnf^ "^^1 ~ = f^Nz = f^N] and the other one with 



hierarchical light and hierarchical heavy neutrinos, and with niy^ ~ eV, rriy^ = \/ '^i^^l, 



This is a reduced version of our more complete work!"! to which we address the reader for 
more details. 

2 Numerical results and conclusions 

For the degenerate heavy neutrinos case, we have gotten rates for all the studied LFV r and fi 
decays that are below the present experimental upper bounds, and we are basically in agreement 
with previous results!^. The largest rates we get, within the explored range of the seesaw 
and mSUGRA parameter space, are for the r decays. Specifically, BR(r fij) ~ 10~® and 
BR(t^ —>■ fi^fi^fi^) ~ 3 X 10~^^, corresponding to the extreme values of tan/3 = 50 and mjy = 
10^^ GeV and for the lowest values of Mq and M1/2 explored. The case of hierarchical heavy 
neutrinos turns out to be much more interesting. First of all, we get much larger branching ratios 
than in the previous case and secondly they are in many cases above the present experimental 
bounds. 
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Figure 1: Dependence of BR(r jj. /i ^■'') with |6l2 1 for arg(S2) = 7r/4, (mjvj , mjvj , mjva ) = (10*, 2 x 10*, 10^"') 
GeV, 61=63= 0, tan/3 = 50, Mo = 400 GeV, M1/2 = 300 GeV, Ao = and sign(^i) > 0. 



For the more interesting hierarchical case we show in Figs. 11151 the numerical results for 
the branching ratios of the LFV r and decays for both types, Ij — s- 3/i and the correlated 
radiative decays Ij — s- li'j. All the LFV r and ^ decay rates are mainly sensitive to tan/3, the 
heaviest neutrino mass rn-TVa, which we have set in these figures to m^Vs = 10^^ GeV, and the 
complex angles in the R matrix 9i and O2, which have been taken in the range 3 < tan/3 < 50, 
< l^il < 3 and < arg(0j) < 7r/4. In Fig. ^ we show separately the various contributions to 
BR(t~ ^~^) as a function of |^2|- The dominant one is the photon-penguin contribution 

(which is undistinguisible from the total in this figure) and the others (Z-penguin, H-penguin 
and boxes) are several orders of magnitude smaller. All the rates for r~ — > /U"/!"//"*" in this plot 
are within the allowed range by the experimental bound, which is placed just at the upper line 
of the rectangle. In contrast, we have obtained that, for the mSUGRA and seesaw parameters 
of Fig. ^ the predicted BR(r —>■ fi'y) is above the experimental bound for the explored range of 
< 1^2 1 < 3. In Fig|2]we show the predictions of BR(/~ — > l~l^lf) and BR(/j —>■ lij) as functions 
of \92\, for all the channels and for different values of arg(02)- It is clear that all the branching 
ratios have a soft behaviour with |02| except for the case of real 62 where appears a narrow dip 
in each plot. In this Fig.|21we see that all the rates obtained are below their experimental upper 
bounds, except for the processes r — > |U7 and /i 67, where the predicted rates for complex 62 
with large |02| are clearly above the allowed region. The most restrictive channel in this case is 
r — > /i7 where compatibility with data occurs just for real 62 and for complex 62 but with |02| 
values near the region of the narrow dip. We also see that the rates for BR{jjL — > 3e) enter in 
confiict with experiment at the upper corner of large |02| and large arg(02) = t^/^- Even more 
interesting are the predictions for BR{1~ —>■ and BR{lj li^y) as functions of \6i\, due 

to the large values of the relevant entries of the Yy coupling matrix. Concretely, can be as 

large as ~ 0.2 for \6i\ ~ 2.5 and arg {61) = 7r/4, and \Y^'^\ and \Yj)'^\ are in the range 0.1 — 1 for all 
studied complex 61 values. The results for BR{lJ —>■ l~l~lf) and BR{lj — > /j7) as functions of 
16*1!, for different values of arg(0i), are illustrated in Fig. 121 We see clearly that the restrictions 
are more severe in this case than in the previous one. In fact, all the rates cross the horizontal 
lines of the experimental bounds except for BR{t~ and BR{t~ — > e~e~e~^). The 

most restrictive channel is now the fi ^ decay. More specifically, we see that all the points 
in the plot of BR{^ — > 67), except for the particular values = and real 61 at the dip, are 
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Figure 2: Dependence of LVF r and fi decays with 1^2 1 with hierarchical heavy neutrinos and complex R, for 
arg(6'2) — 0, 7r/10, 7r/8, 7r/6, 7r/4 in radians (lower to upper lines), {mff-^,mi^2,m]si^) = (10*, 2 x 10*, 10^'*) GeV, 
Oi =03 = 0, tan /J = 50, Mo = 400 GeV, Afya = 300 GeV, Ao = and sign(/i) > 0. The horizontal lines are the 

upper experimental bounds. 



excluded by the experimental upper bound. Also the predictions for BR{^ 3e) are mostly 
excluded, except again for the region close to zero and the dip. We have also explored the 
dependence with the complex ^3 angle, and it turns out that the predictions for all rates are 
nearly constant with this angle. We have also tried another input values for the heavy neutrino 
masses. The results for BR[t — > 3^) are shown in Fig. |1J Here we compare the predictions for 
the three following set of values, {mN^,niM^,mN;) = (10^ 2x 10^ 10^^) GeV, (10^°, 2x 10^°, 10^^) 
GeV and (10^, 2 x 10^, lO^^) QgV. We conclude, that the relevant mass is the heaviest one, mjq.^, 
and the scaling with this mass is approximately as the scaling with the common mass rriTv in 
the degenerate case, namely, {rriNz log Jit-ats)^- Because of this, the rates for the two first sets are 
nearly undistinguisable, and the rates for the third set are about four orders of magnitude below. 
Finally, we consider the very interesting case where the angle ^13 of the Umns is non vanishing. 
It is known that the present neutrino data still allows for small values of this angle, ^13 < 10°. 
The dependence of BR{fi~ — > e~e~e'^) and BR{ii — > 67) with this ^13 is shown in Fig. [^Jwhere 
we explore values in the < ^13 < 10° range. We choose these two channels because they are 
very sensitive to this angle and because they have the most restrictive experimental bounds. For 
this study we assume the simplest choice of i? = 1 , and set the other parameters to the following 
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Figure 3: Dependence of LFV r and /i decays with |6'i| with hierarchical heavy neutrinos and complex R, for 
arg(0i) = 0, 7r/10, 7r/8, 7r/6, 7r/4 in radians (lower to upper lines), (mjvi , m/^j , miVa) = (10*, 2 x 10*, 10^'*) GeV, 
02 = 03 = 0, tan/3 = 50, Mo = 400 GeV, M1/2 = 300 GeV, Ao = and sign(^i) > 0. The horizontal lines are the 

upper experimental bounds. 
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Figure 4: Dependence of LFV r ^ 3/i with |^2| in with hierarchical heavy neutrinos, for different mjv; choices. 
Solid line is for (mAfj , mAfj , mjvg) — (10*, 2 x 10*, 10^*) GeV, dashed line is for (mjvi , mATj , mjva) = (10"''',2 x 
10^°, lO^"*) GeV, and dotted line is for (miVi, mjvj , mATg) = (10*, 2 x 10*, 10^^) GeV. The rest of parameters are 
set to tan/3 = 50, Mo = 200 GeV, M1/2 = 100 GeV, Aq = 0, sign(//) > and arg(6'2) = 7r/4. The horizontal line 

is the experimental bound. 
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Figure 5: Dependence of LFV fi decays with 9i3 in radians with hierarchical heavy neutrinos and R = 1, for 
(mjvi , ?Tijv2 , fTtjVs) ~ (10^,2 X 10*, 10^'*) GeV, tan/3 = 50, = and sign(/i) > 0. The upper lines are for 
Mo = 250 GeV, Ph/2 = 150 GeV and the lower hnes are for Mo = 400 GeV, M1/2 = 300 GeV. The horizontal 

lines are the experimental bounds. 



values: tan/3 = 50, Aq = 0, sign(^) > 0, and {mNi,mN2,rnN3) = (10^,2 x 10^,10^'') GeV. The 
upper lines are for Mq = 250 GeV, M1/2 = 150 GeV and the lower ones for Mq = 400 GeV, 
Ml/2 = 300 GeV. We conclude that, for this choice of parameters, values of 613 in the upper 
region of the interval < ^13 < 10° are clearly disfavored by the data on LFV /z decays. It is a 
quite stricking result. 

In summary, we obtain in the hierachical case much larger rates than in the degenerate 
one, and one must pay attention to these values, because the rates in several channels do get 
in conflict with the experimental bounds. More specifically, the choice of a complex R matrix 
with large modules and/or large arguments of 9i and/or 02 and a light SUSY spectrum is very 
constrained by data. We also confirm that the experimental upper bounds of the processes 
Ij — > kl are more restrictive than the Ij — > l^l^lf ones but all together will allow to extract 
large excluded regions of the mSUGRA and seesaw parameter space. A more precise conclusion 
on the excluded regions of this parameter space deserves a more devoted study. 
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